Boundary Conditions 


The first order approximation for the surface boundary condi- 
tions on pressure and temperature are 


P(m = Mr) = (4.1.3) 


T(m= Mr) =0 (4.1.4) 


Compared to the values for temperature and pressure within the 
star, small values such as P = T = 0 for the surface conditions 
are valid. However, we can improve upon these boundary condi- 
tions by using relations between surface temperature, pressure, 
luminosity, and radius that come from the analysis of stellar at- 
mospheres. 








To estimate the pressure and temperature in the atmosphere of 
a star with luminosity, £7, radius, R, and mass Mr, recall that 
for radiative diffusion 
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where u is the energy density, and we have explicitly labeled 
the frequency dependence of u and «, and £. Since the energy 
density is related to the radiation pressure by 
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(9.1) is equivalent to 
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To simplify this equation, we first get rid of the frequency depen- 
dence by finding a mean value for the opacity. If we write (9.2) 
without the v dependence, 


L Ea dP rad 








Anr2 Rp dr ay) 


then it is clear from a comparison of the two equations that 


k= a Ky Lydv (9.4) 


Note that this is not a Rosseland mean opacity (but often a 
Rosseland mean is used). 


If we now substitute dr = & pdr in (9.3), the equation becomes 
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Compared with the stellar radius, the atmosphere is very thin, 
with r ~ R. Thus we can trivially integrate (9.5) from the stellar 
“surface” (T = 0) to a given optical depth, Tp, by 
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Note here that P.aq(7 = 0) will not be u/3, since that implies an 
isotropic radiation field. (Recall, that in deriving the diffusion 
equation, we assumed that there were just as many particles dif- 
fusing inward across a boundary, as there were particles diffusing 
outward.) If we make the simple assumption that photons only 
diffuse outward at the stellar surface, then P.aq(7 = 0) will be 1/2 
the isotropic value, i.e., P:aa(7 = 0) = u/6. In thermodynamic 
equilibrium, u, = 47 B(T)/c, which means 
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If we now define the effective temperature of an object through 
the equation Lr = 4r R?0T4,, then 7 B(T) = oT 4, and 
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This approximates how the radiation pressure of a star changes 
with optical depth. Similarly, if we relate the radiation pressure 
to the temperature by Paa = aT*/3 (and remember that a = 
4o /c), then 
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Note that (9.8) defines what the optical depth is at the location 
of the photosphere, i.e., at T = Teff, Tp = 2/3. 





Since we now know where the stellar photosphere is in terms of 
Tp, we can use this to compute the pressure at the photosphere. 
In hydrostatic equilibrium, 
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dt (9.9) 





If we integrate this equation from the surface (r = 0) to the 
photosphere (7 = 2/3), and make the (somewhat poor) assump- 
tion that the opacity doesn’t change with optical depth, then the 
photospheric pressure, Pp, is 


E 2GMr 
P 3K R? 





+ P(r =0) (9.10) 


Now let’s assume that at r = 0 the gas pressure is negligible 
compared to the radiation pressure. With this assumption, we 


again have 
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which yields 
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The second term in (9.11) is the ratio of the stellar luminosity 
to the Eddington luminosity. Under most circumstances, Lr < 
4ncGMr/k, so the latter term can be neglected. This leaves 
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Matching Radiative Atmospheres to Interiors 


Although the stellar photosphere boundary conditions represent 
an improvement over the surface boundary conditions, they still 
are not ideal. At the photosphere, the mean free path of a pho- 
ton is long compared to the pressure, density, and temperature 
gradients. This violates the assumptions that went into deriving 
the equations of radiative energy transport. To avoid this prob- 
lem, this atmospheric grid should be used to match the interior 
models below the photosphere. 








To see how this can work, let’s adopt an opacity law of the form 
K = Kop?T? = Ky POT? (9.13) 


where (for an ideal gas) 
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If the envelope is radiative, i.e., 
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which implies 
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This expression explains why one can use surface boundary con- 
ditions (P = 0, T = 0) in a stellar model. If 8 < 4, and a > —1, 


(as is true for most opacities), then as one goes deeper into the 
interior of the star, 
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making P and T independent of the surface conditions. More- 
over, a comparison of (9.14) and (9.16) yields 
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which allows us to relate the temperature to the star’s radius. 
Recall that An GM oT 
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If we assume that the amount of mass contained in surface layers 
is negligible, and use the ideal gas law, then 
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We can thus determine boundary conditions deeper in the star, 
and reach a region where the stellar interior approximations hold. 








or 








Convective Atmospheres 


For cool stars, H~ opacity dominates the surface layers. Since, 
B ~ 9 for this type of opacity, the analysis above does not apply. 
To see what happens, consider the behavior of temperature and 
pressure in a small radiative region on the surface of the star. At 
the surface, the radiative temperature gradient is 
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Now, re-write (9.15) as 
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divide through by T4~, and substitute in Vyaq and Va for P 
and P, 
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The expression for V;aq at some depth in the star is then 
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This can be further simplified if we evaluate V, at the stellar 
photosphere. If we substitute 
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(the latter from (9.12)), and remember that a = 40/c, then 
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This; at =p = 2/3,.and V5. = 1/5: 


Now, consider a cool star, in which H~ opacity in important. At 
low temperatures, x x p!/?T°, and thus a = 1/2 and 8 = 17/2. 
Equation (9.22) for Vraa then gives 
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The equation demonstrates that Vraq will increase rapidly as the 
temperature increases towards the interior of the star. In fact, 
just below the surface, Vraq will become greater than Vaa, Thus, 
the region just below the surface will be convectively unstable. 
This is a feature of all cool stars: if H~ is the dominant source of 
opacity, convection will occur. If we assume that convection will 
begin when Vraa > 0.4, then the temperature at this transition 
zone will be 
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In a similar fashion, we can estimate the pressure at the transition 
zone. If we take (9.21) and substitute for the constant Cı using 
(9.20), then 
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From (9.23) and (9.24), Vp = 1/8 and T/T, = (8/5)?/%, so this 
equation can be directly solved for P. When we plug in the 
numbers for the a = 1/2, 8 = 17/2 case of H~ opacity, then the 
pressure at the transition zone becomes 
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